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Measurements of the vertical velocity of ice crystals us- 
ing a 1.5/im Doppler lidar are described. The statistics 
from a continuous sample of stratiform ice clouds over 17 
months are analysed: the distribution of velocity varies 
strongly with temperature, with average Doppler velocities 
of « 0.2ms"^ at -40°C increasing to « 0.6ms"^ at -10°C 
presumably due to particle growth and broadening of the 
size spectrum. We examine the likely influence of small crys- 
tals less than GOfim by forward modelling their effect on the 
area-weighted fall speed, and comparing the results to the 
lidar observations. The comparison strongly suggests that 
the concentration of these small crystals in most ice clouds is 
much lower than measured in-situ by cloud droplet probes. 
The discrepancy is attributed to shattering of large crystals 
on the probe inlet, and we argue that these numerous small 
particles should not be included in numerical weather and 
climate model parametrizations. 



1. Introduction 

There has been much controversy over the number of 
small sub-60/im particles in ice clouds. Measurements from 
forward-scattering cloud droplet probes suggest that these 
tiny particles are present in very large concentrations (eg. 
Boudala et al 2002, Piatt et al 1989, Gayet et al 2007, Law- 
son et al 2006; see also Heymsfleld and McFarquhar 2002), 
with important implications for the cloud's microphysical 
evolution, its interaction with radiation (Boudala et al 2007, 
De Leon and Haigh 2007, Zender and Kiehl 1994), and re- 
mote sensing (Donovan 2003). On the other hand, there is 
evidence (Gardiner and Hallett 1985, Field et al 2003, Mc- 
Farquhar et al 2007, Heymsfield 2007) which suggests that 
these small particles are in fact instrument artefacts, caused 
by large crystals shattering on the probe inlet, with the de- 
tector then sampling the numerous tiny fragments. 

General circulation models (GCMs) typically represent 
ice particle size spectra using a simple exponential or gamma 
distribution (eg. Wilson and BaUard 1998, MitcheU 1991), 
with the concentration of particles of diameter D given by 



N{D) = iVoI>^exp(-Ai3). 



(1) 
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Figure 1. Example bimodal size spectra for T = — 15°C 
(solid line), — 25°C (dashed line), and — 35°C (dotted 
line). IWC is flxed at O.Olgm^^ 



The parameters A^o and A are either diagnosed as a func- 
tion of temperature T or deduced from the model ice wa- 
ter content (IWC), whilst the parameter ^ is typically as- 
sumed to take a constant value (zero for a pure exponential). 
However, if the measured concentrations of small crystals 
are genuine, equation 1 substantially under-represents their 
numbers. 

To better represent the small crystal observations, 
Ivanova et al (2001) have taken cloud droplet probe and 2D 
optical probe measurements from 17 flights through strat- 
iform ice clouds (a total of over 1000 individual spectra), 
and used them to parametrize a two-mode size distribution 
for input into GCMs. In this scheme, particles larger than 
~ lOO/xm are represented by the simple gamma distribution 
described above, with A diagnosed from the cloud temper- 
ature. However, a tall, narrow bump centred at around 
25/im is also added to represent the small particles seen by 
the cloud droplet probe. Examples of these parametrized 
spectra at temperatures of —15, —25 and — 35°C are shown 
in figure 1 for a fixed IWC of O.Olgm"^. Concentrations of 
crystals 25/im in diameter are typically 2-3 orders of magni- 
tude higher than crystals 250/xm in diameter at these tem- 
peratures, and will strongly influence the lower moments of 
the particle size distribution (eg. optical extinction, depo- 
sition/evaporation rates). A very similar representation of 
the size spectra was suggested for remote sensing studies by 
Donovan (2003). 

Mitchell et al (2008) applied the two-mode parametriza- 
tion outlined above to a GCM and showed that these tiny, 
barely-falling crystals would have a significant effect on ice 
cloud coverage and radiative forcing, particularly in cold 
tropical clouds. It is therefore vital to determine whether 
these concentrations of small crystals are genuine or not, in 
order that simulation of ice cloud in numerical weather and 
climate models is to be realistic. 
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Wcstbrook (2008) suggested that Doppler lidar measure- 
ments could bo a sensitive test for the presence of genuine 
populations of small crystals, since the Doppler velocity 
would be significantly weighted by the small mode (which 
falls at only a few cms~^). In what follows, we test this 
idea, using 17 months of continuous observations from the 
Doppler lidar based at the Chilbolton Observatory in Hamp- 
shire, UK, and compare the results to a forward model of 
Ivanova et al's parametrization of the ice particle size spec- 
tra, both with and without the controversial small mode. 



2. Doppler Lidar observations 

The measurements were made using a 1.5/im Doppler li- 
dar (HALO Photonics Ltd, Malvern, UK). This instrument 
has a maximum range of 10km and measures profiles of 
backscatter and Doppler velocity every 32s, at 36m reso- 
lution, and operated continuously between September 2006 
and January 2008 when the data presented here was col- 
lected. The lidar points directly at vertical: such a con- 
figuration is necessary to measure the fall speeds of the ice 
particles without contamination from the horizontal wind; 
however it can also load to strong specular reflections from 
oriented planar ice crystals, which may bias the Doppler ve- 
locities. This is particularly a problem in mid-level, mixed- 
phase clouds (West brook et al 2009). Comparison of the 
backscatter with a second lidar pointing slightly off zenith 
allows these regions of cloud to be identified and removed 
- see Westbrook et al for more detail. Liquid cloud, pro- 
files with precipitation at the ground, and boundary layer 
aerosol are also removed using the same procedure described 
in that paper. This processing effectively limits our sample 
to stratiform ice clouds which are not raining/snowing at 
the ground, and our statistics should reflect the properties 
of these clouds. A total of 1.1 million 32sx36m pixels were 
analysed; this breaks down into ~ 85000 pixels per 2.5° C 
interval at -10°C, decreasing to « 22000 at -40° C. 

For lidax measurements at non-absorbing wavelengths, an 
ice particle with a given shape and orientation is expected 
to produce a backscatter proportional to its projected area 
A, and the corresponding Doppler velocity would represent 
the area-weighted fall speed of the crystal population. How- 
ever, at 1.5/im there is some absorption as the light reflects 
around the inside of the ice crystal: as the particles be- 
come larger, the amount of absorption increases, reducing 
the backscatter. This leads to the lidar Doppler velocity 
being more strongly weighted towards the smaller particles 
than a simple arearweighting. We can write this as: 



/ N{D)AfvdD 
J N{D)AfdD 



+ w 



(2) 



where w is the vertical wind speed, and the dimensionless 
factor f{D) falls off as a function of crystal size from / « 1 
for small crystals, to / w for very large ones. As a guide, 
a backscattered ray of light with a path length of 200/im 
through the ice crystal gives / = 0.45 (Westbrook et al 
2009). 

Strictly, this analysis is only correct if the particle is much 
larger than the wavelength: for the small crystal mode this 
approximation may not be entirely correct. Mie calculations 
for ice spheres at 1.5//m (not shown for brevity) follow the 
predicted steady roll off in backscatter efficiency due to ab- 
sorption for large particles. For very small spheres there 
is an increasing amount of oscillatory behaviour, leading 
to a slight enhancement in backscatter efficiency at 25/im 
(the peak of the small mode), and a suppressed efficiency 



at 15/im. In practice, for a cloud containing an ensemble 
of different particle shapes, much of this oscillatory struc- 
ture will be smoothed out, and it seems unlikely that these 
departures from the simple geometric optics approximation 
will have a significant effect on (v), or on our conclusions 
below. 

Figure 2 shows the distribution of Doppler velocity for 
all the ice cloud sampled by the lidar over the 17 month 
period as a function of temperature (T is talcen from the 
Met Office 12 km forecast model output over Chilbolton, 
see lUingworth et al 2007). The distribution is broad, and 
has a clear temperature trend with particles falling faster 
at warmer temperatures, indicating the influence of particle 
growth/aggregation. Note the sign convention of negative 
Doppler velocities for particles falling toward the lidar. 

Smoothing the lidar data over 10 minute periods before 
calculating the statistics leads to a narrower distribution as 
shown in the right hand panel of figure 2, removing small 
scale variations in w, although some influence from lower fre- 
quency waves may still be present in the distribution. Since 
we are looking at stratiform clouds we expect that large 
scale ascents will be weak. Westbrook (2008) estimated that 
the maximum fall speed for a compact 60/nm ice particle is 
w 0.06ms~^ , so given the observed velocity data we conclude 
that the small crystal mode does not completely dominate 
the lidar signal in most clouds; however there could be a 
more subtle influence. We investigate this by comparing 
the observed mean values to the forward modelled velocities 
from Ivanova et al's two- mode size spectra. 

3. Comparison with ice parametrizations 

3.1. Including small crystals 

We have calculated the mean Doppler velocity for all ice 
clouds sampled as a function of temperature from the distri- 
bution above, and this is plotted in figure 3. As noted above, 
the temperature trend is marked, with velocities increasing 
from « 0.2ms"^ at -40°C to « 0.6ms-i at -10°C. For 
this average over 17 months of data, we expect any residual 
influence of w in the average values to be positive and very 
small compared to the total observed Doppler velocities in 
figure 3. 

Ice particle size spectra were calculated as a function of 
temperature as described in Ivanova et al (2001), superim- 
posing two gamma distributions: 

N{D) = NosD''' exp(-A.D) + NotD"' exp(-A^I>). (3) 

For the large particles fii = 2.64 and is diagnosed from T 
using fits to the aircraft data. For the small mode, Hs = 3.24 
and As is fixed at 3.65 x 10'*m~'^ (giving a mode centred at 
approximately 25/im). There is a simple partitioning of IWC 
between the two modes: this depends weakly on Xe, and was 
inferred from aircraft data. For the range of temperatures 
sampled here, the fraction of the total IWC contained in the 
small mode was essentially constant at 11%. This percent- 
age fixes the ratio Nos/Not, and the concentration of small 
crystals relative to the large mode. The spectra in figure 1 
are examples of this prescription. 

These size spectra were then used with the mass- and 
area-diameter relationships assumed in that paper to calcu- 
late the area-weighted fall speed of the crystals using the 
formula in Mitchell (1996). Although Westbrook (2008) 
shows that Mitchell's formula can overestimate the veloc- 
ity of small particles, switching to his formula for the small 
mode had a minimal effect in this case. Air densities were 
assumed for each temperature beised on a standard mid- 
latitude atmosphere (International Organization for Stan- 
dardization 1975). Because A (rather than iVo) is diagnosed 
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Figure 2. Filled contour plot showing distribution of lidar Doppler velocity. Left panel shows distribution of raw 32-s data; 
Right panel shows same when 10 minute smoothing is applied. Grey scale indicates frequency in a particular velocity and 
temperature bin; frequencies are normalised such that the total in each temperature bin = 1. Bins are 0.05ms~^ x 2.5° C. 
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Figure 3. Average lidar Doppler velocity as a function of temperature (black line with markers). Left panel shows 
comparison to forward modelled velocities from Ivanova et al including the small mode: solid grey line is prediction when 
the particles are assumed to be planar polycrystals; dashed grey line shows the sensitivity of changing this assumption 
to aggregates; dot-dash shows change to rosettes. Right panel shows forward modelling of large mode only (black dash). 
Results of using A(T) relationship measured by Piatt (1997) are shown by dot-dash curve; grey shading shows prediction 
from Wilson and Ballard's scheme for IWCs in the range 0.001 — O.lgm"^. 



from T, the forward modelled velocity for a given tempera^ 
ture docs not depend on IWC. 

The forward modelled velocity is plotted as a function of 
temperature in figure 3: given the analysis in section 2 we 
expect this to be faster than the measured lidar velocities 
due to the absorption effect at 1.5/jm. However for tempera- 
tures warmer than w — 35°C we find the opposite trend: for- 
ward modelled velocities are slower than the observations, 
by almost a factor of two at warm temperatures. We also 
note that the forward modelled curve has a much weaker 
temperature trend than the observations, varying from ap- 



proximately 0.3ms~^ at -40°C to 0.35ms"^ at -5°C. These 
discrepancies both suggest the small crystal mode is having 
too strong an influence on the forward-modelled velocity. 

One concern about the above analysis is that it may hinge 
on the assumption of a particular set of (perhaps unrepre- 
sentative) mass- and area-diameter relationships. In their 
parametrization, Ivanova et al assume the particles to be 
planar polycrystals. However they also consider the influ- 
ence of changing this assumption to bullet-rosettes (leading 
to a slightly altered partitioning of IWC). Using the rela- 
tionships in Mitchell (1996) we have forward modelled the 
result of this alternative assumption, and this is plotted in 
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figure 3. Now the forward modelling leads to even slower ve- 
locities, in the region 0.2rns^"'; again the temperature trend 
is much weaker than observed. Finally, it could be argued 
that in many clouds the large particle mode is typically dom- 
inated by aggregates: we have investigated the effect of this, 
again using the relationships given by Mitchell (1996), and 
assuming the same partitioning of IWC as for planar poly- 
crystals. The forward modelled fall speeds are very similar 
to the planar-polycrystal curve. Since analysis of aircraft 
images (Korolev et al 2000) indicates that most ice particles 
in stratiform clouds are irregular polycrystals or aggregates 
these three curves should give a good idea of the sensitiv- 
ity to realistic changes in particle shape, and so we can be 
reasonably confident that the inconsistency of the forward 
model with the observations is due to the representation 
of the particle size spectrum, and not the representation of 
particle shape. 

Given figure 3 it seems clear that the small mode ex- 
erts too rrmch infiuence on the parametrized size spectra for 
temperatures warmer than ~ — 35°C, and this adds weight 
to the suggestion that the measured concentrations of small 
crystals are likely to be artefacts most of the time. To inves- 
tigate this further, we have also forward modelled the size 
spectra with the small mode removed. 

3.2. Large mode only- 
Figure 3 shows results of forward modelling Ivanova et 
al's size spectra if only the large mode is included (dashed 
line). The velocities are faster, and at warm temperatures 
are close to the Doppler lidar observations, whilst falling 
somewhat faster than the lidar velocity at cold tempera- 
tures. Also shown is the result of altering the A(r) relation- 
ship to that measured by Piatt (1997) in midlatitude ice 
cloud: this leads to somewhat faster fall speeds, lying to the 
left of the lidar observations at all temperatures; the slope 
of the curve with temperature is similar to the observations. 

For comparison we have also forward modelled velocities 
from the parametrization of Wilson and Ballaxd (1998). In 
this scheme a simple exponential size distribution (equation 
1 with fi = 0) is assumed, with No diagnosed from tempera- 
ture and A deduced from the model IWC. Mass and velocity- 
diameter relationships are provided in Wilson and Ballard, 
but no A{D) relationship is defined, so we apply the rela- 
tionship for aggregates/polycrystals given in Mitchell (1996) 
to calculate the area-weighting. The grey shaded region in 
figure 3 shows the range of velocities predicted for IWCs be- 
tween 0. 001-0. lgm~^, and this velocity range is consistent 
with the lidar observations. The temperature trend also ap- 
pears to be realistic given that higher IWCs are correlated 
with warmer temperatures. 

At present it is not possible to evaluate which of these 
parametrizations is best: if it becomes possible to parame- 
terise f{D) in a sufficiently robust manner, we could forward 
model the true lidar Doppler velocity (rather than the area 
weighted velocity), and such a comparison would place a 
stronger bound on which relationship is most realistic. 

4. Discussion 

Doppler lidar observations have been used to infer 
whether measurements of high concentrations of small crys- 
tals in ice clouds are likely to be genuine or not. Frequency 
distributions of measured Doppler velocity show a system- 
atic increase in fall speed with temperature, due to broaden- 
ing of the size distribution. Forward modelling of the area- 
weighted fall speed from two-mode size-spectra is not consis- 
tent with the lidar data and strongly suggest that the large 



rmmbers of small ice particles measured iii-situ are likely to 
be instrument artefacts in many cases. Forward modelling 
with the large mode only gave results which were much more 
consistent with the observations. We therefore argue that 
inclusion of a small mode in numerical weather and climate 
model parametrizations is likely to be unrealistic. However 
we are cautious not to polarise the issue: clearly in a number 
of situations small ice particles are present, especially near 
cloud top where there may well be recently frozen droplets 
or haze particles present, and crystals in the early stages 
of growth. It would also be valuable to make similar mea- 
surements in cold tropical clouds where one might expect 
high concentration of crystals to be nucleated. However, we 
conclude that the measurements of numerous small parti- 
cles at all heights in mid-latitude ice cloud are inconsistent 
with our lidar observations, which are of course unaffected 
by shattering. 

One limitation of this study is that the lidar is unable 
to penetrate optically thick cloud: this means that in deep 
systems, only the lowest 2km or so will be observed. This 
may mean that at cold temperatures the data is biased to- 
wards thinner clouds producing less ice (which the lidar can 
see through), and presumably lower fall speeds. This would 
imply that the average lidar velocity at cold temperatures 
may in fact be even faster than in figure 3 and this would act 
to further reinforce our conclusion about the small crystal 
mode. Another limitation of this work is that 25% of the ice 
cloud sampled was rejected because of specular reflection 
from oriented planar crystals, predominantly in mid- level 
clouds. One might worry that this could bias the sample; 
on the other hand, co-incident radar reflectivity measure- 
ments show a similar distribution for both specular and 
non-specular cloud (West brook et al 2009). More investi- 
gation of the potential influence of this removed cloud on 
our statistics would be valuable. 

Future work will focus on calculation of f{D) using ray 
tracing sinmlations of polycrystal and aggregate shapes; this 
will then allow us to completely forward model the true li- 
dar velocities from model parametrizations. We also note 
that implicit in equation 2 is the assumption that the parti- 
cle shape (or ensemble of different shapes) does not change 
systematically as a function of particle size. Since we have 
taken a large sample of different clouds we might expect 
any such correlation to average out; theoretical sensitivity 
studies would be valuable to back this up. 

For models which diagnose iVo as a function of T and pre- 
dict A from IWC, coincident radar reflectivity (Z) measure- 
ments may allow the parametrizations to be tested, since Z 
is related to A through A^o . The combination of sinmltaneous 
radar reflectivity, radar Doppler spectra and lidar Doppler 
velocity has the potential both to test model parametriza- 
tions as a whole, and also to help discriminate which individ- 
ual elements of that parametrization are realistic and which 
are not. We aim to persue this approach in future work to 
evaluate the Met Office Unified Model cloud scheme. 
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